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Abstract The actual research deals with determining 
by a new protocol the necessary parameters consider- 
ing a three-dimensional model to simulate in a realistic 
way the turning process on machine tool. This paper is 
dedicated to the experimental displacements analysis of 
the block tool / block workpiece with self-excited vibra- 
tions. In connexion with turning process, the self-excited 
vibrations domain is obtained starting from spectra of 
two accelerometers. One three axes accelerometer placed 
on the tool and one unidirectional accelerometer placed 
on the front of bearing of spindle. The existence of a 
displacements plane attached to the tool edge point is 
revealed. This plane proves to be inclined compared to 
the machines tool axes. This plane contains an ellipse 
that is the place of the points of the tool tip displace- 
ments. We establish that the tool tip point describes an 
ellipse. This ellipse is very small and can be considered 
as a small straight line segment for the stable cutting 
process (without vibrations) . In unstable mode (with vi- 
brations) the ellipse of displacements is really more vis- 
ible. A difference in phase occurs between the tool tip 
displacements on the radial direction and on the cutting 
one. The feed motion direction and the cutting one are 
almost in phase. The values of the long and small ellipse 
axes (and their ratio) shows that these sizes are increas- 
ing with the feed rate value. A weak growth (6%) of the 



long and small axes ratio is obtained when the feed rate 
value decreases. The axis that goes through the stiffness 
center and the tool tip represents the maximum stiffness 
direction. The maximum (resp. minimum) stiffness axis 
of the tool is perpendicular to the large (resp. small) el- 
lipse displacements axis. The self-excited vibrations ap- 
pearance is strongly influenced by the system stiffness 
values, their ration and their direction. FFT analysis of 
the accelerometers signals allows to reach several im- 
portant parameters and establish coherent correlations 
between tool tip displacements and the static - elastic 
characteristics of the machine tool components tested 
(see 0|). 
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Depth of cut (mm) 

Large (small) axis of ellipse attached with 
the cutting tool tip displacements (m) 
Large (small) axis of ellipse attached with 
the forces charge points (m) 
Block Tool 
Block Workpiece 
Feed rate (mm/rev) 
Number of revolutions (rpm) 
Vector attached to the components system 
stiffness 

Normal direction of the plane P u 

Unitary vector, support of the ellipse axis i 

(i = a, b) situated in the plane P u 

Plane attached with the cutting tool edge 

displacements 

Cutting edge radius (mm) 

Cutting tool edge radius (mm) 

Total roughness (fim) 
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Nomenclature 



T Time that corresponding to one workpiece 

revolution 

u Tool edge point displacement (m) 

x (z) Radial (axial) direction 

y Cutting feed axis 

a Clearance angle (degree) 

a K,(xy) Angle of the main stiffness direction in the 
plane (x,y) (degree) 

a K,(yz) Angle of the main stiffness direction in the 
plane (y,z) (degree) 

At Temporal dephasing between two signals 

ip u Dephasing between the displacements com- 
ponents attached to the tool edge point (de- 
gree) 

7 Cutting angle (degree) 

A s Inclination angle of edge (degree) 

K r Direct angle (degree) 

O e (xy) Main displacements direction angle of the 
tool edge point in the plane (x,y) 

Q e (yz) Main displacements direction angle of the 
tool edge point in the plane (y,z) 



1 Introduction 

The machine tools dynamic phenomena come from the 
elastic system machine / cutting process interaction |17j , 
[19] , [2"T] . [2"5] . This interaction is thus the dynamic as- 
pects generating source classically met in the machine 
tools use [18], [25] • The cutting process actions applied 
to the elastic machining system cause relative tool / 
workpiece displacements that generate vibrations [20] . 
These modify the chip section, the contact pressure, the 
speed of relative movement etc. Consequently, the cut- 
ting process instability can cause the dynamic system 
instability of the machine tool: vibrations appear [34] . 
They have undesirable effect on the machined workpiece 
surface quality and on the tool wear [2], [37]. Mainte- 
nance problems may occur, even ruptures of the machine 
tools elements [35] ■ Thus, it is necessary to develop nu- 
merical and experimental models ensuring to study the 
vibratory phenomena during machining, in order to ob- 
tain the stable process in cutting [15] , [21] , [28] , [43] . 

Considerable efforts were developed to model the cut- 
ting process correctly. However, at the present time the 
suggested solutions are still far from providing sufficiently 
relevant and general models to translate in an accept- 
able way all the available experimental results [10] , [H] . 
Models become very complex for studying the three- 
dimensional model and nonstationary cutting conditions, 
particularly in the case of the vibratory situations, [3J, 
[T4] . [18] . [31], [33], [H]. A model for the expression of 
cutting forces uncertainty in different cutting conditions 
was proposed [1] . The model enables the uncertainty pre- 
diction of cutting force measurements, for a defined cut- 
ting parameters range, based on calibration data and 



measurements collected in a few experiments run using 
a factorial design. The results showed that the contri- 
butions in the uncertainty budget from single cutting 
force measurements are balanced between the calibra- 
tion error on the channel where the force is measured, 
and the cumulative error due to the cutting parameters 
variation. Being responsible for approximately 50% of 
the uncertainty in single force experiments, the cutting 
parameters variation may be responsible for the scatter 
of cutting force measurements found in the literature. 

Today, the machine tools have rather high dynamic 
behaviours with regard to rigidities and damping capac- 
ities. However, the causes of the vibrations appearance, 
or unstable modes, are given by the the cutting dynam- 
ics process under various working conditions. In other 
words, under certain conditions, the cutting interaction 
process with the machine tool elastic system causes the 
appearance of vibrations. Self-sustained vibrations are 
observed. These vibrations, due to the cutting process, 
have frequencies close to the natural system frequencies, 
and are generated by the cutting actions variation, vari- 
ation which depends on the various parameters charac- 
teristic of the process [12] , [13] , [20] . 

These parameter modifications lead to the variations 
of relative displacements of tool / workpiece and tool / 
chip. Variations of the chip section are generated, induc- 
ing variations of the cutting actions, and thus maintain- 
ing the vibrations. It is thus advisable to take into ac- 
count these parameter variations in a three-dimensional 
model cutting process which requires a specific deter- 
mination. These problems were already studied in the 
past. Different methods were used to estimate workpiece 
clastic deflections under cutting forces in turning by us- 
ing analytical modelling of workpiece elastic deflections, 
by finite element (via Castigliano's second theorem) [2], 
and by finite-differences [38]. In [28], a real-time com- 
pensation scheme is developed based on a cutting force 
induced error model. 

In this paper, an experimental study is developed to 
determine and characterize a series of parameters neces- 
sary for the dynamic model of cutting process in agree- 
ment with static aspects obtained previously [4]. After 
having exposed the tests protocol used (see Sect. [2]), dis- 
placements are measured and a frequency analysis is car- 
ried out (see Sect. [3]). A displacement plane attached of 
the tool point (see Sect. [4| is highlighted starting from 
the accelerometer data. The deep analysis included in 
Section [5] shows that these displacements describe an el- 
lipse. The evolution of the ellipse axis ratio can be char- 
acterized by a slightly increasing function depending on 
the feed motion while their ratio follows an opposite law. 
This ellipse is located in a tilted precise plane compared 
to the machine tool axes. Before concluding, Section [6] 
gives a correlation between the tool tip stiffness and dis- 
placements on the one hand, and on the other hand be- 
tween stiffness center [3] and the central axis of the dy- 
namic cutting process. 
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2 Test protocol 

2.1 Approach 

In turning operation, the workpiece has a rotary motion 
whereas the tool has a linear translation. It is obvious 
that both the tool and workpiece vibrate in material re- 
moval. Although workpiece vibrations impact both cut- 
ting instability and product quality including surface fin- 
ishing, most models developed for investigating surface 
roughness [55], [ID], [H] do not consider workpiece vi- 
brations at all. 

A new tests protocol is here proposed to determine a 
series of led phenomcnologic parameters during the cut- 
ting process in particular in the self-sustained vibrations 
case. According to [4], the machining system static be- 
havior presents couplings in all the directions. It is thus 
advisable to expect that it is the same in dynamics. 

The approach adopted for this study constitutes a 
part of a global vision of research presented in the Fig. [TJ 
As in the static case, the dynamic characterization is car- 
ried out starting from the interaction of two elements: 
block tool (BT) and block workpiece (BW) forming the 
machining system and knowing that, a priori, each el- 
ement taken separately have a different vibratory be- 
haviour (cf. figure [3]), [5]. 
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Fig. 1 General approach considered 



by piezoelectric sensors adapted for faster acquisitions 
and being adequate and compatible with the vibratory 
phenomena necessary to measure. 

The three-dimensional dynamic character is highligh- 
ted by seeking the various existing correlations between 
the various parts of the machining system and the vari- 
ous parameters evolutions, which ensure to characterize 
the process. 




Cross slide rest 



Six component 
dynamometer 



Fig. 2 Experimental machining system and metrology envi- 
ronment in turning process 



2.3 Data acquisition protocol 

The self-excited vibrations characterization method, is 
based on various work [26], [35], [39] and on the results 
of the preliminary dynamic study undertaken in [5], is 
summarized in Fig. [3] 



2.2 Experimental device 

Dynamic cutting tests are carried out on a conventional 
lathe (Ernault HN 400) for which the spindle speed docs 
not exceed 3 500 rpm. The main components of the test 
machining system used are presented without its data- 
processing environment in the Fig. [21 The machining 
system behavior is analyzed through one three-direction 
accelerometer fixed on the tool and using two unidi- 
rectional accelerometers positioned on the lathe, on the 
front bearing of the spindle, to identify the influence of 
this one on the cutting process. Moreover, a six-components 
dynamometer [16], being used as tool- holder [46], is po- 
sitioned on the lathe to measure all the cutting forces 
(forces and torques). This study is presented in the pa- 
per [5]. The dynamometer used is an evolution of |16j . 
The active part made of the strain gauges is replaced 



Different components of 
machine tool system 



Dynamo meter 



>BT 



Fiiang plate 



Carnage 



cutting process: two kinematics chains are running: spindle chain and advance movement 



dnve motor and kinematics chain of the spindle are turning 



No cutting process: only dnve motor is turning 



Fig. 3 Natural frequencies domain of different components 
of the machining system 
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The static study undertaken here (cf. |4] and [6]) 
allows to identify perfectly the machine tool frequency 
domain and each one of its components under various 
operating conditions. Moreover, a virgin frequency band 
at any natural frequency was not found. The machine 
tool was turned on without part to machine. Thus, the 
appearance of frequencies in this zone when the cutting 
process exists, highlights the vibratory phenomenon in- 
duced by this one. In what follows, the cutting param- 
eters are selected constant, excluding the depth of cut 
(ap). 

On the test machining system, the instantaneous spe- 
ed of rotation is controlled by a rotary encoder directly 
related to the workpiecc. The connection is carried out 
by a rigid steel wire, which allows a better transmission 
of the behavior (Fig. [2|) . During the cutting process, the 
number of revolutions is controlled constantly at nearly 
690 rpm and a negligible variation cutting speed of about 
1 % is detected. 

All measurements are collected on a PC and are thus 
treated in real time. For these tests, the used tool was a 
noncoatcd carbide tool (TNMA 16 04 12) without chip 
break geometry. The cutting material is a chrome molyb- 
denum alloy type (AiSi 4140). The cylindrical test-tubes 
have a diameter of 120 mm and a length of 30 mm 
(Fig. [4|). The dimensions retained for these test-tubes 
were selected using the finite elements method coupled 
to an optimization method using SAMCEF® software 
presented in [BJ. 
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Table 1 The tool geometrical characteristics 




Fig. 5 The tool geometrical details. 



The machining system is presented as an elastic system 
and the interface BW / BT is the closure of this system 
(Fig. [BJ. The vibrations appearance is strongly condi- 
tioned by the elastic structure behaviour of the studied 
system [7|. 



Fig. 4 Tested workpiece 



-- (gj Woik-piece | 

I Tool * J 



BLOCK TOOL x 



BLOCK WORK-PIECE 



Moreover, the geometry of the tool (Fig. [5]) [27] is 
characterized by the cutting angle 7, the clearance angle 
a, the edge angle of inclination X s , the direct angle K r , 
the nozzle radius r e and the sharpness radius R. The tool 
insert is examined after each test and is changed if nec- 
essary of wear along the cutting face (Vb < 0.2 mm ISO 
3685), that may disturb the studied phenomenon. The 
tool characteristics used are presented in the Table [TJ 

Remember that the assembly of BW detailed in [BJ is 
considered as a block of interdependent elements, which 
is subjected to the vibrations due to the cutting process. 



mtim u '////////// 

Fig. 6 Model of machining system and BW - BT interaction 



Next research step is to pass to the precise analysis 
of the different parameters measured. 
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3 Tests and analyses 

3.1 Experimental conditions 

Many tests are carried out under the conditions specified 
in the previous paragraph. They arc considered for the 
depth of cut ap = 1 mm, 2 mm, 3 mm and 5 mm in order 
to compare signals and to dissociate the machine vibra- 
tion influence of those generated by the cutting process 
itself. 



3.2 Workpiece/Tool/Chip displacements 

The signals treatment of the three-direction accelcromc- 
ters gives, by integration, displacements according to the 
three directions. An example is presented in the Fig. [7] 
related to a cutting depth of 2 mm. For these tested 
conditions, the system is stable; the signal has very low 
amplitudes, on the order of micrometer; and the chip sec- 
tion remains constant (Fig. [7] and [8]) . The part surface 
quality is correct, with a total roughness (i? t ) of 4.3 fim. 



1-10* 

3- 10-= 

6-10- 5 

4- 10-s 
- 2-10- s 
| 

-2-10-5 

5 

-4-10- s 
-6-10- s 
-8-10-5 




Fig. 8 Chip related to ap = 
spindle speed N = 690 rpm 



3 mm, f = 0.1 mm/rev and the 




Fig. 9 Chip related to ap = 4 mm, 
the spindle speed N = 690 rpm. 



f = 0.1 mm/rev , and 



TiiiK [s] 



Fig. 7 Displacements components signals of the tool tip for 
a turning operation: ap = 2 mm, f = 0.1 mm/rev, and the 
spindle speed N = 690 rpm 



The analysis of measurements for the depth of 3 mm 
cut, give similar results; the displacements amplitudes 
are in the same region, lower than 10 micrometers, and 
the chip section does not vary (Fig. [H]) ; this one has the 
same continuity characteristics as for ap = 2 mm (Fig. [7]), 
but very light undulations appear on the workpiece sur- 
face. 

Actually, the characteristics start to change distinctly 
from the depth of 4 mm cut. The tool point displace- 
ments amplitudes increase, and the signal sinusoidal char- 
acter has high amplitude compared to the previous cases. 
The system starts to become unstable. The chip develops 
important periodic undulations and significant section 



variations appear (Fig. [9]). All chips are type 1.3 (ISO 
3685). The machined surface part has weak undulations. 

The increase of cutting depth until 5 mm allows to 
reach a mode of totally unstable cutting and the self- 
excited vibration appearance is clearly identified for ap 
= 5 mm fFig.llip. Around this working point, many tests 
are carried out in order to determine the self-excited vi- 
brations extent. These tests are carried out maintaining 
the same ap value and the same workpiece revolutions 
number N but for various feed motion f, illustrated in 
the Table [2] 



ap(mm) 


N (rpm) 


f(mm/rcv) 


5 


690 


0.05 | 0.0625 | 0.075 | 0.1 



Table 2 The cutting process parameters 



It is important to specify that, for the four feed val- 
ues used (Table [2]), the cutting process was in unstable 
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mode, more or less marked, with chip section variations 
and very different surface undulations on the workpiccc 
for each imposed feed value. In order to reach the best 
possible precision, we then study the extreme case ob- 
tained with the maximum values which involve the most 
important dynamic effects: cutting depth ap = 5 mm, 
and feed motion value f = 0,1 mm/rev, but maintaining 
a constant cutting speed. 

In Fig. [10] are presented the three signals related to 
the tool tip displacements components according to the 
three directions of workspace. 

% 



MO- 4 i "r 




Zoom figure 1 1 



Fig. 10 The tool tip displacements on the three directions 
x, y, and z related to the machine tool and considering ap = 
5 mm, f = 0,1 mm/rev, N = 690 rpm 

A signals zoom related to the tool tip displacement 
components is presented in the Fig. [JTJ as well as the 
results of cutting process vibrations on the workpiccc 
and chip. 

In this situation, the chip presents important tight- 
ened periodic undulations. The section variations seem 
to be in the case of a discontinuous cutting. The work- 
piece surface quality becomes poor and the undulations 
left on the workpiece after the tool passage are very im- 
portant (Fig. [TTj) . The chip characterizations and the 
workpiccc are detailed in [3]; in particular, the undula- 
tions and the chip section variations are quantified there. 



3.3 Displacements frequency analysis related to the tool 
tip 

The tool tip displacements components frequency anal- 
ysis considering the three directions related to the ma- 
chine tool is carried out starting from a Fast time Fourier 
Transformation (FFT). An FFT spectrum example of 
the three acceleration components delivered by the three- 
directional accelerometer placed on the tool is given in 
the Fig. QT] The frequencies above 1 000 Hz are elimi- 
nated knowing that, according to the literature, the self- 
excited vibration frequency domain is situated between 
120 Hz and 200 Hz [23], [2B], [3jJ, [32], [35], [42]. 




Fig. 11 Signals zoom related to the tool tip displacements 
components according with cutting parameters: ap = 5 mm, 
f = 0,1 mm/rev., N = 690 rpm 



The FFT analysis of several tests shows that the ac- 
celeration frequencies of vibrations during the cutting 
process are about 150 Hz to 200 Hz. They are slightly 
lower than the machine natural frequencies (Fig. [3]). How- 
ever, we know that the self-excited vibrations appear 
near the natural and the machining system are not de- 
pending on the number of workpiccc revolutions [24j . 

On the example of vibrations measurements during 
the machining, presented in the Fig. 1121 it should be 
noted that the cutting vibrations fundamental frequen- 
cies are around 190 Hz for the three axes, with a domi- 
nating amplitude on the y axis (cutting axis, according 
to Fig. [5]). These frequencies are lower than those ob- 
tained ([5]) during the machining system dynamic char- 
acterization with a no-cutting process, which are located 
at approximately 550 Hz (Fig. [3]) [5] . The recorded fre- 
quencies, therefore, correspond well to the self-excited 
vibrations related to the cutting process. The frequency 
peak around 190 Hz is found for the other tests, and it 
is constant according to the feed motion in turning. 

Another very important parameter for the self-cxcitcd 
vibrations is their amplitude. It depends on the energy 
quantity introduced into the system. Therefore, the am- 
plitude is influenced on the one hand by the system elas- 
tic parameters, and, on the one other hand, by the cut- 
ting parameters: cutting speed, feed motion and depth 
of cut. The parameter that presents the most important 
influence is the depth of cut. The vibrations amplitude 
increases in a quasi linear way according to the increase 
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Fig. 12 FFT spectrum corresponding of the accelerations 
acquires on the three directions; cutting parameters selected 
ap = 5 mm, f = 0.1 mm/rev. and N = 690 rpm 



depth of cut until it exceeds the elastic system stabil- 
ity limit in agreement with literature data published by 
Kudinov l26l. 



The determination of this plane is presented in ap- 
pendix (Sect: 18. 2p . The existence of this plane is con- 
firmed by the displacements signals observation fFig. fTTj) . 
It is henceforward a question of locating this plane in 
three dimension space. 



4.2 The tool tip displacements plane localization 

The tool tip displacements plane P u is characterized by 
its normal vector noted n u (Fig. [14]) . 



4 Tool point displacements plane 



4.1 Acceleration analysis of vibration data 

After integration, the acceleration data analysis allows 
to establish the existence of an attached displacements 
plane of the tool tip, which, describes an ellipse (Tig. [F3")) . 
Certain parameters are a slightly increasing depending 
of the feed motion, as we will demonstrate further. 



Fig. 14 The tool tip displacements plane P u : projection of 
the normal vector n u on the axis related to the machine tool 




Fig. 13 The tool tip displacements plane P u 



The characterization of this plane by its normal is 
carried out for each test depending on the feed motion 
in turning (Table [3]). 



f(mm/rev) 


0.05 | 0.0625 | 0.075 | 0.1 


Normal 


n u 


x direction 


-0.071 


-0.071 


-0.058 


-0.056 


y direction 


-0.186 


-0.186 


-0.206 


-0.216 


z direction 


0.98 


0.98 


0.977 


0.975 



Table 3 The normal n u of the plane P u attached of the tool 
tip displacements on three directions of cutting, depending 
on feed motion 



We note that the normal component according to the 
feed motion direction almost does not evolve according 
to this one. On the other hand, a very light decrease 
(respectively growth) according to the radial direction 
(respectively, the axis of cutting) is observed in depen- 
dence with the feed motion values. These variations are, 
however, so weak that they do not affect the standard 
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of the plane direct normal which thus can be consid- 
ered constant according to the feed motion. The exis- 
tence and the determination of this plane provide sim- 
plifying important informations about the configuration 
to be adopted in order to write a semianalytical model 
that is under development [4], [II]. In particular, we 
will further see that this plane allows to bring back the 
cutting three-dimensional problem, with three directions 
displacements, using a simpler plane model that is in- 
clined compared to axes (x, y, z) related to the machine 
tool. 



5 Displacements analysis of the tool point 

5.1 Stable and unstable process comparison 

In the tool tip displacements plane P u , it is possible to 
establish that the tool tip point describes a very small 
ellipse in the stable case cutting process (without vibra- 
tions); this ellipse could be assimilate with a small seg- 
ment of straight lines (ap = 2 mm), according to the 
Fig. 1151 a. On the other hand, in unstable mode (with vi- 
brations), the ellipse (according of tool tip displace- 
ments (ap = 5 mm) is really more visible (Fig. 1151 b). 



= 2 mm 





Fig. 15 Comparison between the stable (a) and unstable 
cutting process (b) 



5.2 Phase control between the tool tip displacements 
components 

A phase difference is occurring between the tool tip dis- 
placements following the x radial direction and the y 
cutting direction or the z feed motion direction (these 
last two being almost in phase). This phase difference 
is evaluated starting from the comparison of two signals 
that have different phases. More precisely, phase differ- 
ence <p u is evaluated starting directly to the temporal 
signals and using the relation: 



where T is the time corresponding to a phase differ- 
ence of 360° and At to temporal phase difference between 
the two compared components signals. Fig. [16] presents 
an example of phase difference <p u between radial dis- 
placement x and the two displacements components fol- 
lowing y and z; they are almost in phase. A phase dif- 
ference of 28° is thus obtained for the studied case, ap 
= 5 mm and f = 0.1 mm/rev. 




At 

±2tt— , 



(1) 



Fig. 16 Example of phase difference between x radial tool 
tip displacements and other displacement components (fol- 
lowing cutting direction and feed motion direction). 



During the experimental tests according to different 
feed motion values, the variation of this phase difference 
is very weak. For the value f= 0,05 mm/rev, the phase 
difference remains equal to the phase difference obtained 
for the feed value f = 0,1 mm/rev. 



5.3 Tool tip displacements ellipse approximation 

The ellipse is characterized in detail in the appendix (cf. 
Sect: I8.3P on the plane P u , with the normal direction 
n u , by its large axis a u and its small axis b u determined 
using the method of least squares, in the reference system 
(n ua , n u b) unit vector related to the ellipse axes (Fig. 

HZ). 

The values of a u , b u and their ratio a u /b u are pre- 
sented in the Table [H The value analysis of the large 
and small ellipse axis shows that these two sizes are in- 
creasing with the feed value. On the other hand, a weak 
growth (6%) of the ratio a u /b u is obtained when the feed 
value decreases. The ellipse geometrical proportions are 
thus quasiconstant according to the feed weak evolution 
chosen here. 



6 Discussion on static and dynamic aspects 

6.1 Stiffness / displacements correlation 

Taking into account the machining system static study 
realized in [6], it is possible to check that the stiffness 
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Fig. 17 Coordinate system changing to obtain the ellipse 
attached to the tool point displacements 



f (mm/rev) 


a u (m) 


b u (m) 


a u /b u 


0.1 


7.71. 10~ b 


1.72. lO" 5 


4.48 


0.075 


5.6.10" 5 


1.23.10" 5 


4.59 


0.0625 


3.87.10" 5 


0.83.10" 5 


4.66 


0.05 


2.83.10" b 


0.51. 10" 5 


4.76 



Table 4 Values depending on the feed for large and small 
ellipse axes that correspond to tool tip displacements for ap 
= 5 mm, f = 0.1 mm/rev 



principal directions angles are comparable with the el- 
lipses angles attached to the tool tip displacements (Ta- 
ble ©. 



static case 


dynamic case 


a K{yz) 






"e(xy) 


f (mm/rev) 






77" 


69 U 


0.1 


76° 


65° 


77° 


69° 


0.075 






79° 


67° 


0.0625 






79° 


65° 


0.05 



Table 5 Comparison values between the principal directions 
angles of the tool tip displacements and the stiffness principal 
directions of the machining system 



Between the vector n f (characterizing the stiffness 
principal directions) and the normal vector n u situated 
in the attached tool tip displacements plane, another cor- 
relation exists. Indeed, calculation shows that the angle 
between these two vectors does not exceed 2.5°, in our 
case (Table [6]). It is advisable, moreover, to note that the 
difTcrcncc between the two vectors is slightly increasing 
depending on the feed value. 



f(mm/rev) 


arccos ( .. "f "" ... ) 

\ ||nrll-l|n u ||| J 


0.1 


2.5" 


0.075 


1.9" 


0.0625 


1" 


0.05 


0.8" 



Table 6 Correlation between nf - stiffness principal direc- 
tions and n u - normal direction of the plane P u attached of 
the tool tip displacements 



6.2 Correlation between tool displacements / stiffness 
center 

By comparison with the static analysis results developed 
in [4], it is highlighted a coupling between the system 
clastic characteristics BW and the vibrations generated 
by cutting process. As we could expect it, the self-excited 
vibrations appearance is strongly influenced by the sys- 
tem stiffness values, their ratio and their direction. 



Axis ofstifliiess centre passing by tool point 




x (m) 

Fig. 18 Correlation between stiffness center and the ellipse 
attached of the tool tip displacements 



The stiffness center (CR) is given in [5] (and a sum- 
mary of the stiffness center determination is given in ap- 
pendix l8TTj) . The axis, which passes by the stiffness center 
and the tool tip, represents the maximum stiffness direc- 
tion. This maximum stiffness axis is drawn in Fig. 1181 It 
is orthogonal with the normal of the plane P u , and it is 
also perpendicular to the large ellipse axis. This explains 
well why the greatest displacement takes place, logically, 
on the axis where the stiffness is minimal. This correla- 
tion validates the experimental protocol established [5] 
for this case of turning. 
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To continue, the P u displacements plane determina- 
tion and the elliptic characterization of this plane motion 
provide significant information on the configuration to be 
adopted writing a simplified three dimensional model. It 
is clear at the end of this research that the system must 
be considered according to this plane axes and not ac- 
cording to those of the machine tool. 



6.3 Correlation between stiffness center and central axis 
of the dynamic process in turning 

As example, for a bearing test on 68 revolutions (ap = 
5 mm, f = 0.1 mm/rev. and N= 690 rpm) we calcu- 
lated the beam of the central axis of the tip tool point 
small displacements torsor, their intersection point and 
the stiffness center corresponding to this configuration. 
All these results are presented on Fig. [19] The line that 
binds the stiffness center (CR) and the point of central 
axes intersection passes by the tool tip point. 



Place of the central 

axes 




Fig. 19 Central axis / stiffness center 



tool tip displacements place show the coherence of the 
presented new approach and consolidates the experimen- 
tal procedures implemented. 

From the present work, the following conclusion can 
be drawn: 

- The acceleration data analysis allows to establish the 
existence of an attached plane of the tool tip displace- 
ments tilted compared to axes related to the machine 
tool. 

- In this plane, the tool tip displacements describes an 
ellipse. 

- We have established that the tool tip point describes 
a very small ellipse (assimilated with a straight-line 
small segment) in the case of the stable cutting pro- 
cess (without vibrations); in unstable mode (with vi- 
brations) the ellipse of displacements is really more 
visible. 

- A phase difference is occurring between the tool tip 
displacements following the radial and cutting direc- 
tion or the feed motion direction (these two last being 
almost the same phase). 

- The long and small ellipse axis values show that these 
two sizes are increasing with the feed value. 

- A weak growth (6%) of the ratio long on small axes 
is obtained when the feed value decreases. The ellipse 
geometrical proportions are, thus, quasiconstant ac- 
cording to the feed weak evolution chosen here. 

- The axis which passes by the stiffness center and the 
tool tip represents the maximum stiffness direction. 

- The maximum (resp. minimum) stiffness tool axis is 
perpendicular to the large (resp. small) ellipse dis- 
placements axis. 

- The self-excited vibrations appearance is strongly in- 
fluenced by the system stiffness values, their ratio and 
their direction. 

It is also important to note for this research the plane 
determination attached to the tool tip displacements which 
is carried out according to an ellipse the surface of which 
is an increasing function depending on the feed motion 
value; this aspect is in perfect coherence with the power 
level injected in the system during machining. 



7 Conclusion 

Compared to other authors who use only the tool stiff- 
ness [T2|, [IS], [35], this study takes into account, more- 
over, stiffness resulting from the system fixture/tool- holder 
(BT). The important information on the static stiffness, 
obtained in [6J, certainly makes it possible to determine 
the privileged displacements directions. 

The stiffness center and the rotation center were ob- 
tained experimentally [6J . The direction of minimum dis- 
placement was also defined starting from the experimen- 
tal model. The confrontation of these results: stiffness 
center, principal stiffness directions, central axis and the 



8 Appendix 

8.1 Reminder of stiffness center determination 

It is well known that any machine tool is characterized 
by principal directions of deformation. These principal 
directions of displacement are a function of the machine 
structure, its geometrical configurations, and cutting pa- 
rameters used. We can observe either a very stiff behav- 
ior, or very rubber band according to the direction. Here, 
the tool is regarded as forming integral part of block 
tool BT. The aim is to determine the stiffness center 
CRbt of elastic system BT [30]. This stiffness center 
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corresponds to the block BT rotation center compared 
to the lathe bed. Obtaining CRbt consists in finding 
the intersection points of different perpendicular to dis- 
placements. 

The stiffness center is obtaining based on the tool 
point O that is origin of the coordinate system (x, y, 
z). The procedure to obtain the stiffness center is based 
on the virtual work principles and detailed in the figure 
l20l the imposed charge follows all three directions, Fi 
being known; (i = x, y, z). Measurement of the displace- 
ment vector (linear and rotation displacement) under the 
charge di,j is obtained in two points Pi,j in each direction 
(i = x, y, and z) and (j = 1 and 2). The direction (Dij) 
is released being straight line containing dij vector that 
includes the point Pi,j . 



Fy 




Fig. 20 Experimental procedure to obtain stiffness center 
CRbt 



Using the method of least squares we can minimize 
the distance dx between the lines. Each of points Mi on 
the lines (Dij) are separated by a minimum distance \ii 
= 1,8 mm and maximum deviation Si = 2°. 

In the next step, three mean plans Pi are denned by 
approximation of the points Mj and containing the lines 
Dij. We can draw the normal rii of each plan Pi that 
include the point Mi. 

We obtain a linear system with three vectors' equa- 
tions and three unknown factors, after that is seeking to 
obtain the intersection point CRbt using least squares 
minimization method. In reality, the directions of these 
three normal lines cross (nearly at r = 1.2 mm) in only 
one point noted CRbt which corresponds with stiffness 
center. 



8.2 Determination of the attached plane of the tool tip 
displacements 

The modeling of the attached plane of the tool tip dis- 
placements starting from the experimental results is car- 
ried out using the software Mathcad©. We seek to de- 



termine the plane that passes by the points cloud repre- 
senting tool point displacements (Fig. n~3|) : 

ax + by + cz + d = 0. 

Using the parameter e rr for the errors, the following 
is obtained : 

e rr (x, y, z, Xp : yp : Zp : x n , yn: z n ) — 
[M(x, y, z) - P{x pi y pi z p ).n u (x n , y ni z n )} , 

where: 

M = {x, y, zf , P {x p , y p , z p }* , n u = {x n , y n , z n } , (2) 

in this equation, superior index t indicates the trans- 
position operation. Let us express the E rr function us- 
ing e rr and introduce the components of displacements 
(u Xi ,u yi ,u Zi ) on the three space directions, it comes: 

E rr (x, y : Z : Xp, 7/p, Zp : X n7 yn: Z n ) — 

N 

^ ^ &rr {^Xi j ^yt : j 2Cp : Up 3 Zp , X n , Vn : Zn): 
i=0 

and : 

EN t-^ n X^ N 

i=0 Ux * ., _ l~,i=0 U Vi _ 1^,1=0 U Zi 

Xp ~ N lVp ~ N ' p ~ N 

The vector V is computed now by minimization: 

V 771271 (E rr: Xp 7 yp 7 Zp 7 X n: Vn: Zn) , 

where: 



V = 



r vi 






= Vp: 


« v * 


Zp: 




— %n : 


v 5 


= Vn, 


{ v 6 


Z n . 



It results from it that the n u normal components in 
the plane of the tool tip displacements are given by: 

{ Va, 
n u = { V 5 , 

(v 6 . 

For the present case study in the section lOl (Tablel3j) . 
with ap = 5 mm and f = 0.1 mm, it comes: 

( -0.056, 
n u = { -0.216, 
[ 0.975, 
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the equation of the searched plane is then : 

P u {s,t) = V p + S.Ui +t.u 2 , 

where s and t are constants, and V p is equal with: 



V P = < v 2 , 



Vi, 

v 3 , 



Ul = 



and Ui, U2 are given by: 
u A n u 



-,u 2 = nAu L 



u 



n A u u | 

Uq being the plane orientation vector: 



1, 
1. 
1. 



i if 





I 



X (in) 



"I - 

■Mtr* 



-1 10'' 



-0 



--i lct 



Fig. 21 2D projected ellipse attached of the tool tip point 
displacements 



8.3 Ellipse approximation 

Using the ellipse plane determination, it is possible to 
build the matrix M rep expressed by: 



ui! u 2l n u i 

[M rep ] = Ui 2 U 22 n ul 

Ui 2 u 23 n ul 
Another matrix is defined by 



[M U ] = [U}-[P], 

with [P] defined by the equation (2) and [U], being 
the vector of the tool tip point displacement given by: 



[M r 



cos(9 a ) — s'm(8 a ) 
sin(0 Q ) cos(9 a ) 
1 



The angle 6 a is obtained using minimization method 
and Mathcad© software. We obtain: 



[M„ 



[M ro 



\M„ 



and after that is obtained the tool tip displacement 
related to the base system (cki, aa): 

[U 2 ] = [M U3 ]. 

Consequently, the ellipse equation is: 



P] 



(u ai ) 2 (u a2 ) 2 
a 2 b 2 

using parameters f = 0.1 mm/rev and ap = 5 mm. 



Now, we determine the new ellipse configuration cen- 
tered on the origin of the basic reference system (x, y, Acknowledgements The authors would like to thank the 
z): Centre National de la Recherche Scientifique (UMR 5469) for 

the financial support to accomplish the project. 

[M U2 ] = [M rep ] 4 .[M u ]. 



The new displacements of the tool tip are : 



[Ui] = [M U2 ]. 

By changing the coordinate system, we project the el- 
lipse in the new two-dimensional reference system fFig.[2T 
(011,0(2) using the rotation matrix [M rot ] 
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